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Abstract: We study the properties of the holographic CFT dual to Gauss-Bonnet gravity in 
general D (> 5) dimensions. We establish the AdS/CFT dictionary and in particular relate 
the couplings of the gravitational theory to the universal couplings arising in correlators of 
the stress tensor of the dual CFT. This allows us to examine constraints on the gravitational 
couplings by demanding consistency of the CFT. In particular, one can demand positive en- 
ergy fluxes in scattering processes or the causal propagation of fluctuations. We also examine 
the holographic hydrodynamics, commenting on the shear viscosity as well as the relaxation 
time. The latter allows us to consider causality constraints arising from the second-order 
truncated theory of hydrodynamics. 



Contents 



|]. Introduction 

||]. Gauss-Bonnet gravity 

|. AdS/CFT dictionary 
|3.1| Central charge Ct 



3.2 Holographic calculation of t% and £4 

3.3 Constraints 



4.1 



Gauss-Bonnet black hole perturbations 

Scalar channel 
Shear channel 
Sound channel 



4.2 



4.3 



Causality constraints 
^TT| Shear channel 
Sound channel 
Scalar channel 
Plasma instabilities 



5.2 



5.3 



5.4 



Second-order hydrodynamics 
^Tl| Holographic Gauss-Bonnet hydrodynamics 
6 . 1 . 1| Speed of sound 



6 . 1 . 2| Shear viscosity-to-entropy ratio 



6.1.3 Relaxation time and causality violation 



|7|. Discussion 

[X]. Energy fluxes in terms of the three-point couplings 

[B]. Conformal tensor fields 

Q, Useful expressions 

[D|. Coefficients for the sound channel 




I 

12 




n 

16 

n 

19 

2C 

n 

22 

25 

1 

2C 

27 
28 
35 
39 
41 
42 



- 1 - 



1. Introduction 



The AdS/CFT correspondence [[[], ^ [| or more generally gauge/gravity dualities provide 
a theoretical framework in which to study (certain) strongly coupled gauge theories. For 
example, this approach allows for the calculation of transport coefficients for gauge theory 
plasmas at strong coupling by means of relatively simple supergravity computations, while 
these calculations are prohibitively complicated by any other conventional methods — for 
example, see @, [|, gj. Recently this topic has been of great interest motivated by the discovery 
of the strongly coupled quark-gluon plasma (sQGP). In particular, it was found that the ratio 
of shear viscosity to density entropy of any fluid dual to Einstein gravity is precisely 1/Air 
[0, H 1 . Despite the fact that these holographic calculations deal with gauge theories which 
are quite exotic compared to QCD, this result still seems to come remarkably close to the 
value measured for the sQGP |To[ ]. 

Originally it was conjectured that these holographic calculations provided a universal 
bound: rj/s > I/An |?J. However, it is now accepted that this conjectured bound is violated 
in string theory by the effect of higher curvature interactions in the gravitational action. Still 
the precise string theory constructions where these higher curvature terms are under control 



only allow for perturbative violations of the bound [11, 12, 13 1. It is certainly also of interest 



to explore situations where finite violations of the bound occur. A useful framework for such 



explorations was found to be Gauss-Bonnet (GB) gravity [14, 15]. The original studies were 
made with five-dimensional GB gravity theory which is dual to a four-dimensional CFT but 
this analysis can easily be extended to any D > 5. Some work in this direction already 
appears in |l(| |l7| . The present paper provides a comprehensive study of holographic GB 
gravity in arbitrary dimensions. 

An overview of the paper is as follows: We begin with a brief review of Gauss-Bonnet 
(GB) gravity coupled to a negative cosmological constant in section ||. In section ||, we 
investigate the AdS/CFT dictionary for these gravitational theories in an arbitrary number 
of dimensions. In particular, we calculate the central charge Ct appearing in the two-point 
function of the stress tensor and the parameters £2 and t\ appearing in the energy one-point 
function describing certain scattering experiments, first proposed in [^] . With these results, 
we determine the constraints on the GB coupling arising from the requirement that the energy 
flux in the experiments is everywhere positive. In section ||, we construct the equations of 
motion for metric perturbations propagating in a black hole background. These equations are 
then examined in section || to study causality violations in the dual CFT. We find that the 
constraints imposed on the GB coupling to avoid such acausality precisely match the positive 
energy flux constraints derived in section [3|. We examine holographic hydrodynamics for GB 
gravity in section |(| In particular, by studying the propagation of sound waves in the dual 
plasma, we derive the relaxation time, as well as the ratio of the shear viscosity to entropy 
density. Here we also consider causality constraints within the framework of second-order 



1 Recently this universality has been extended to T — which is described by extremal black holes in the 
bulk §. 
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hydrodynamics. We conclude with a brief discussion of our results in section ^. We also have 
some appendices containing various technical details. In particular, appendix [A| describes the 
calculation relating the scattering parameters ti and £4 to the couplings A, B and C which 
determine the three-point function of the stress-energy tensor. 

While we were in the final stages of preparing this paper, ref. []l9| appeared which also 
explores causality constraints in Gauss-Bonnet gravity. 



2. Gauss-Bonnet gravity 

Consider GB gravity in D > 5 spacetime dimensions, defined by the following action 

(D- !)(£>- 2) 



1 

w 1 ^ 



d D x, 



-9 



where 



A4 — Rabbit 



L 2 



abed 



+ R + 



L 2 A 



GB 



(Z>-3)(£>-4) 



4R ab R ab + R 2 



(2.1) 



(2.2) 



Of course, this curvature-squared interaction is precisely the Euler density of four-dimensional 
manifolds and so it does not effect the gravitational equations of motion unless D > 5. The 
solutions describing planar AdS black holes take the form p0| 



ds 1 



D-2 



I? dr 2 



i=l 



r 2 f(r) ' 



(2.3) 



where f(r) is given by 



f±(r) 



2A 



1 ± 



1 -4A f 



„D-1 



(2.4) 



In fact, in the following, we will only consider the solutions / = /_ (r) as these will be the 
only ones to correspond to nonsingular black holes in a ghost-free vacuum |21, ^2|. Note that 
in this class of solutions, the horizon appears at r = r+. Using the definition 



foo = lim f(r) 



1 - VI - 4A C 



(2.5) 



' ~ - v 2 A, ; :• 

we have normalized the coordinates above so that lim^oo gu/gxx = — 1- This choice was 
made to set the speed of light to one in the boundary metric {i.e., in the dual CFT). Further, 
by setting r+ = 0, we recover the AdS vacuum metric in Poincare coordinates. Examining 
g rr , we can see that the AdS curvature scale L is related to the parameter L in the action as 
L 2 = L 2 1 foo- In the following, we restrict our discussion to A GB < 1/4 since no AdS vacua 
exist for larger values of A GB as can be seen, e.g., from eq. (|2.5|). 

The Hawking temperature of this black hole solution is given by 

r+ D — l 



^L 2 vTo 



(2.6) 
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The energy and entropy densities are simply calculated as 

D-2 r?- 1 _ D-2 /4vr f^,\ D ~ 2 , , ...„_, 

(2.7) 
(2.8) 

Further note that we find that e oc T d and s oc T^ 1 as expected for a CFT in d = D — 1 
dimensions (in the absence of a chemical potential). Further, these expressions satisfy the 
precise relation e = ^^Ts, again as expected for a conformal plasma. 

3. AdS/CFT dictionary 

In this section, we develop the dictionary relating the couplings in GB gravity theory ( |2.1| ) to 
parameters which characterize the dual CFT. Since we are only dealing with the gravitational 
sector of the AdS theory, we are looking to examine the behaviour of the stress energy tensor 
of the CFT. So for example, the central charges, c and a, appearing in the trace anomaly 
are universal parameters characterizing a four-dimensional CFT p3| and can be calculated 



in holographic context [24]. For GB gravity, these calculations were performed for d = 4 in 
[25 and for d = 6 in [17]. However, while the trace anomaly calculations can be extended to 
examine CFT's in higher dimensions, the details change in each dimension and the number 
of spacetime dimensions must be even. Hence we do not pursue this approach here. 

However, there is a "central charge" common to CFT's in any number of dimensions, d. 
This is the coefficient characterizing the leading singularity in the two-point function of two 
stress tensors [26, 27] : 2 

( T ab {x) T cd (0) > = ^ X abM {x) , (3.1) 

where 

l a b,cd( x ) = 2 CW^O-TwOc) + I a d{x)hc{x)) - -iVabVcd , (3.2) 

and 

Iab(x) = T]ab ~ 2^^ . (3.3) 

This structure is completely dictated by the constraints imposed by conformal symmetry 
and energy conservation p6| , f27|| . Of course, in four dimensions, this coefficient is related 
to the standard central charge c which appears as the coefficient of the (Weyl) 2 term in the 
trace anomaly: Ct = (40/7T 4 ) c. Hence Ct will be one of the coefficients which we calculate 
holographically in the following to establish our AdS/CFT dictionary for GB gravity in general 
dimensions. 



2 Here and throughout the following, we are assuming a Minkowski signature for the metric. Note that in 
eq. (3.3), x a = rj ab x b (i.e., x = -t). 
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To extend the dictionary further we need to identify additional parameters that play an 
analogous universal role for CFT's in any number of spacetime dimensions. With our focus 
on the stress energy tensor, the obvious next step is to look for universal parameters in the 
three-point function, as was extensively studied for CFT's in general dimensions by [26. 27]. 



There it was shown that conformal symmetry was powerful enough to determine the form 
of the three-point function up to five constants, which are labeled A, B, C, V and £ in @. 



Conservation of the stress-energy imposes further constraints which allow us to reduce the 
number of independent parameters to three with: 3 

£ = (d ^ i)A+ 'Mpl B . d Ji±Mc. (3.4) 

Further one finds that Ward identities relate the two- and three-point functions and so Ct 
can be expressed in terms of the parameters characterizing the three-point function |26| , p7[ : 

_ Q d -i (d-l)(d + 2)A-2B-A(d+l)C , . 

° T ~ ~2 dJdT2) ' (3 ' 5) 

where &d-l = 2-K d l 2 /T(d/2) is the area of a unit (d-l)-sphere. In fact, one can perform a 
holographic calculation of the three-point function f28|| , however, extending these calculations 
to GB gravity would require an exhaustive and exhausting analysis. Therefore, we choose an 
indirect route to determining these coefficients in the following. 

In particular, we will consider extending the analysis of the energy flux or "energy one- 
point functions" in |J| for CFT's in an arbitrary spacetime dimension d. This approach is to 



consider an "experiment" in which the energy flux was measured at null infinity after a local 
disturbance was created by the insertion of the stress tensor eij T tJ . The energy flux escaping 
at null infinity in the direction indicated by the unit vector n then takes the form 



d-2 



(3.6) 



where E is the total energy. The structure of this expression is completely dictated by the 
symmetry of the construction. Hence two coefficients, t% and £4, are constant parameters 
that characterize the underlying CFT. Note that the (negative) constants appearing in the 



3 They also find that this general framework meets exceptions in low dimensions, with two independent 
parameters in d = 3 and one, in d — 2. These reductions arise because various tensor structures that are 
independent for d > 4 are not independent with a small number of dimensions. 

4 One comment is that this general discussion only applies for d > 4. It should perhaps be evident that 
d — 3 is an exception since in this case, there are not enough spatial directions to consider rotations in the 



space orthogonal to n. In fact, we find that the starting point (3.6) is not quite correct since 



for d = 3 . (3.7) 
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two factors multiplied by £2 and £4 were chosen so that these factors contribute zero net flux 
when integrated over all directions. The negative sign of these constants leads to interesting 



constraints on the coefficients £2 and £4, which we discuss below in section 3.3. 

As presented in pjj, £2 and £4 can straightforwardly be determined by a holographic 
calculation and we generalize these calculations to GB gravity in any number of dimensions 



in section 3.2 below. As discussed in [18], the energy flux is directly related to the three-point 
function and hence the two coefficients £2 and £4 can be determined in terms of A, B and C. 
A lengthy calculation presented in appendix |A| yields 

_ 2(d+l) (d-2)(d + 2)(d + l)A + 3d 2 B-4d(2d + l)C 
2 ~ d (d-l)(d + 2)A-2B-A(d+l)C ' ( j 

(d+1) (d + 2)(2d 2 -3d-3)A + 2d 2 {d + 2)B-4d(d+l)(d + 2)C 
U ~ d (d - l)(d + 2) A - 2 B - 4(d + 1) C ' 



As further discussed in |18[, a nonvanishing £4 in a four-dimensional CFT implies the 
action of the dual gravity theory must contain terms cubic in the Riemann tensor. This 
analysis readily extends to any number of dimensions and so since such interactions do not 
appear in the GB gravity action ( |2.ip studied here, the holographic CFT must have £4 = 0. 
Combined with eq. (|3 . 9[) , the vanishing of £4 imposes the constraint 



(d + 2) (2d 2 - 3d - 3) A + 2 d 2 (d + 2) B - 4d (d + 1) (d + 2) C = (3.10) 
for the theories studied here. 
3.1 Central charge Ct 



First, we perform a holographic calculation of the central charge Ct appearing in eq. (3.1). 



We follow closely the derivation of the two-point function given in 2£j. We consider metric 
fluctuations propagating in the AdS^+i vacuum geometry. It is convenient to write the latter 
as ds 2 = L 2 /u 2 (rj ab dx a dx b + du 2 ) where as above L 2 = L 2 //^. 5 We choose a gauge where the 
perturbation components 5g uu = 5g au = at the AdS boundary. If we write the remaining 
components as 5g a b = L 2 /u 2 H a b, then the (on-shell) quadratic action for H a f, reduces to the 
following boundary term 

fd-1 r 
h = r^zr(l " 2A GB ioo) / d d xu l - d H ab d u H ab , (3.11) 

8£p JdM 

where the indices are simply contracted with 5 a b- Imposing the boundary conditions 

H ab (u = 0,x) = H ab ( X ), (3.12) 



Hence the most general expression for d = 3 is simply 

E 



(£(n)) d=3 = — 



l + U 



eijn l n 3 \ z 1 



(3.8) 



3 Note that u — (LL)/r where r is the radial coordinate used in section 
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the full bulk solution for can be written as 

Hab(U ^ ] = J=k / ^ + |x - x'|Y lab ^ ~ X ' } ' (3 - 13) 



where T a b,cd precisely matches the tensor structure appearing in eq. (3.2). Note that Sg uu and 
5g au are also nonvanishing in the bulk |^9| but we ignore these polarizations because they do 
not contribute in the quadratic action ( |3.11| ). The quadratic action now becomes 

h - (1 - ^ d _ 1 Jdxdy . (3.14) 

Varying the above expression with respect to Hat, then yields the two-point function of the 
dual stress tensor and upon comparing with ( |3.1[) , we find 



d+1 r[d + i] 

° T ~ d-l^T[d/2] ^ f ] ' ( } 



3.2 Holographic calculation of £2 and £4 

In this section, we perform a holographic computation of the energy flux (|3.6| ). We follow 



closely the approach presented in [25| and so only sketch the salient steps of the calculation. 



We are interested in determining the energy flux 

(O|Ot g ( n )0|Q) 

{£{n)) = (ojotojo) • (3 - 16) 

In the present case, the state is being created by a tensor insertion with O ~ eijT lJ (x) where 
the polarization carries only spatial indices and eij(x) oc e~ lEt . Hence we see that this flux 
is determined by the three- and two-point functions of the stress tensor in the CFT. This 
expression will depend on the polarization tensor e« and the unit vector n l indicating the 
direction in which the flux is measured. As noted above, the SO(d — 1) invariance of the 
constructions fixes the final result to take the form given in eq. ( |3.6| ). 

If we adopt coordinates x^ = x° db x d ~ l , the energy flux measured at future null infinity 
in direction specified by n r can be written as 

£(n) = Jim^J °° dx~(x + ) 2 (l + (t^(x + , x~, n) + Tl j (x + , x~, n)j nP . (3.17) 

We have introduced index notation where i = 1, . . . , d — 2 while i,j = 1, . . . , d — 1 as usual. 
The superscript x on the stress tensor simply indicates that these operators are defined in flat 
space with coordinates x a . This notation is useful because next we introduce new coordinates 
y a following @: 

V + = -^T, y-=x--^, y» = -^. (3.18) 
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Note that the original Minkowski metric described by x a is conformal to a flat space metric 
described by y a : 

ds 2 = -dxUx- + (dx*) 2 = ~ dy+d f + {dyi? ■ (3-19) 



The utility of transforming to y a is that null infinity, x + — > oo, is now mapped to the 
(hyper)plane y + = 0. In this plane, we also have 

V 1 = i-r • (3-20) 



In terms of the new y coordinates, the energy flux (3.17) becomes 



" + 00 

£(n) = n d - 1 I dy-T y __(y + = 0,y-,y l ) (3.21) 



oc 



with Q = 1/(1 + n d x ). We have also introduced the y superscript on the stress tensor to 
indicate the conformally transformed operator: 

^=^%%^- < 3 - 22 > 

To compute (£ (n)) holographically, we must turn on graviton perturbations whose bound- 
ary values source the appropriate energy-momentum tensor insertions. As the first step 
[18, |3(|, we consider the following shockwave background 6 



T 2 / \ T 2 dn 2 

ds 2 = -( 6(y+)W(u, y*)(dy + ) 2 ~ dy + dy- + (dy*) 2 + —— 2 . (3.23) 

In solving for the shockwave profile, the full nonlinear equations of motion for GB gravity 
reduce to a simple linear equation for W(u, y l ) and in fact the curvature-squared terms do 
not contribute |3(], . The solution is chosen with the appropriate asymptotic behaviour to 
source the operator £{n) ~ J dy~ T : 



W(u, j/*) = N W Q d - 1 -r-r (3.24) 



d/2 

~d~ J 

u + (y l — y' 1 ) 2 

with y H = n l /(l + n d ~ l ) as above and Nw is a normalisation constant. 

It remains to add metric perturbations corresponding to the operators O in eq. fl3.16Q . 
For simplicity, we choose a particular polarization with e\2 = a and all other components van- 
ishing. This description applies in the x coordinate system and the dual metric perturbation 
5g a b = L 2 /uH a b would have the boundary condition 

H x i x 2(u = 0) = ae~ iEt = a e'^ (3.25) 



s In this case, u = L 4 /r 2 where r is the radial coordinate introduced in section]^. 
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and all other components vanishing. Next we wish to transform to the y coordinates and 
examine the overlap of the graviton with the shockwave at y + = 0. Changing the coordinates, 
it can be shown O, |25| that the relevant graviton profile is given by 



Hyi y 2(y + =0,y ,y\u) 



1 



(3.26) 



In general, one also includes additional graviton polarizations to ensure the perturbation is 
transverse and traceless in the bulk, however, they will not contribute to the three-point 
function in GB gravity [25]. 

To find the three-point function, we add this perturbation to the metric ( |3.23|) and 
evaluate the action on-shell. We are led to examine terms proportional to W(H y i y 2 ) 2 . After 
integration by parts and using the equations of motion, the relevant part of the cubic effective 
action reduces to the following boundary integral: 



ft 1 



d d y du H yW d 2 _H yly2 W(u, y l ) <%+) 1 



+ 



' 2/ooAgb 

2Aqb/oo 

(D-3)(D-4) 



In the above expression, 



d\W + d%W -4d u W 



W 



T 2 ). (3.27) 



(3.28) 



After inserting the solution ( 3.24 ) for W(u,y z ), we obtain 

T 2 = 2(D-l)(D-2) 



D 



(3.29) 



To normalize our result, we must divide by the two-point function ( T\ 2 T\ 2 ), which was cal- 
culated in the previous section. Upon fixing the normalisation Nw of W(u, y l ) appropriately, 
we find that the result takes the desired form (13.61) with 



i+t 2 



n\ + n 2 



1 



d- 1 



where the coefficient t 2 is given by 



t 2 



4/ooA G b (D-1)(D -2) 



1-2/ooAgb (D-3)(D-4) 
Implicitly, as expected, we have also found that = for GB gravity. 

3.3 Constraints 



(3.30) 



(3.31) 



Because of the negative constants appearing in eq. ( |3.6| ), it is easy to see that if the coeffi- 
cients t 2 and t4 become too large, the energy flux measured in various directions will become 
negative. Avoiding this problem then imposes various constraints, as in 18]. To extract the 
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various distinct constraints, we first fix the unit vector, e.g., n l = (1, 0, 0, • • • ). We then orga- 
nize the polarization tensors e^- according to their rotational properties under the SO(d — 2) 
group that leaves n l invariant. There are three possibilities and each produces a distinct 
constraint: 



• Tensor (spin 2), e.g., £23 = £32 = a and all other components vanish, 

1 2 



1 



t 2 



u > 



d - 1 d 2 - 1 

• Vector (spin 1), e.g., E\ 2 = £21 = a and all other components vanish, 

1 



1 



d-l 



t-2 



d 2 



2 \ 1 

— tA+-t 2 >o 



Scalar (spin 0), e.g., £ij = a x diag(— (d — 2), 1, 1, • • • ), 



d-l 



t-2 



d 2 - 1 



d 



d-l 



(t 2 + U) > 



(3.32) 



(3.33) 



(3.34) 



These results provide a simple extension to a general dimension d of the constraints for four- 
dimensional CFT's derived in 

Using the expressions in eq. ( |3.9| ), the above constraints ( |3.32| ), (|3.33| ) and ( |3.34| ) can be 
translated to constraints on the parameters A, B and C, e.g., , using eqs. ( |A.17|) and ( |A.18| ): 



Tensor: (d - 2) (d + 2) A + 2 dB - AdC < , 
Vector : (d - 2) (d + 2) A + (3d - 2) B - 8 d C > , 
Scalar : i3 - 2 C < . 



(3.35) 
(3.36) 
(3.37) 



Note that the three-point couplings obey an additional constraint arising from the uni- 
tarity of the CFT. The latter implies that Ct is a positive quantity and so from eq. ( |3.5| ), we 
have 

(d- l)(d + 2)A-2B-A(d+l)C > 0. (3.38) 



This inequality was already assumed in deriving the previous constraints (3.35-3.37). 

Before turning to the implications for the dual GB gravity, one might consider the results 

17 



for free fields [27, 26 



A 
B 
C 



1 



n 3 



o 3 



d 3 



'/is 



(d-l)3 

(d - 2)d 3 



d 3 
d-3 

d 2 



(d-i) 3 n * + Y n f + 



(d-2)d 3 
d-3 



(3.39) 



2j2 



(d-2) 2 d 



A(d-l) 3 ' Us+ 4 n/+ 2(d-3) 



(d - 2)d 3 „ ' 



7 The general results for the tensor fields are derived in Appendix [b| 
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where Qd-i is the area of a unit (d-l)-sphere, as defined below eq. (|3.5|) . Also n s indicates 
the number of real conformal scalars and hf is the number of (massless) fermionic degrees of 
freedom. Hence hf = 2 L rf / 2 J f or a massless Dirac fermion in d dimensions — or in the notation 
of [26|, hf = tr(l) where tr is the Dirac trace. Finally for the case that d = 2n, we have 
also included the possible contribution of (n-l)-form potentials, for which the standard free 
kinetic term is also conformally invariant. Then ht denotes the number of degrees of freedom 
contributed by these tensors. Generally, for a single (n-l)-form potential in d = 2n, we would 
have ht = T(2n — 1)/T(n) 2 . For example, with d = 4, this would just be an Abelian vector 
field with two degrees of freedom, i.e., ht = 2. For n odd, we might also consider constraining 
the tensor by demanding that the field strength be (anti-)self-dual, in which case the previous 
result for h t would be multiplied by 1/2. 

At this point, it is interesting to evaluate the constraints ( 3.35| - ^37| ) above with these 



free field results (|3.39| ) for A, B and C 



id 2 - A)d 3 

Tensor: (d - 2) (d + 2) A + 2 dB - 4 dC = — ^—h t , (3.40) 

Vector: (d - 2) (d + 2) A + (3d - 2) B - 8 dC = -{d + 2)d 2 hf , (3.41) 

Scalar: B - 2C = - ^ ~ 4 -*f n s . (3.42) 
2(a — iy 

Hence with various combinations of the free fields, we are able to precisely fill out the allowed 
region that is defined by requiring a positive energy flux in eq. (3. 6). 8 For example, with 



ht = 0, we reach the boundary defined by the tensor constraint in eq. ( |3.35| ), i.e., (d— 2) (d + 
2) A + 2dB — AdC = 0. This boundary surface is mapped out by allowing n s and hf to vary 
but note that eq. ( p. 41 ) shows that hf = corresponds precisely to the line where this tensor 



boundary intersects that set by the vector constraint (|3.36 ). Similarly, from eq. (|3.42 ), we 



see that n s = corresponds to the intersection of the boundaries set by the tensor (|3.35| ) and 
scalar (|3.3?1 ) constraints. 

Now using the AdS/CFT dictionary established above for GB gravity, we may translate 
the constraints above to constraints on the GB coupling A GB . We begin with the constraint 
Ct > and comparing eq. fl3.15| ), we find 

1 - 2/ooA G b > . (3.44) 



However, given the definition of /oo in eq. (2.5), we find that this constraint is always satisfied - 



assuming A GB < 1/4, otherwise no AdS vacua exist. This stems from choosing only to consider 

8 Momentarily, we are treating n s , hf and h t as continuous variables here. This may be seen as a convenient 
approximation in the regime where Ct 2> 1. As a further note, recall that as discussed in footnote ^, the flux 



(3.8) only contains a term in d = 3. In this case, there are not enough spatial dimensions to establish the 



tensor constraint (3.32). The vector and scalar constraints reduce to 

-4<t 4 <4. (3.43) 



With d = 3, the free field expressions (3.39) yield £4 = 4(n s — hf)/(n s + hf) and hence the free theories again 
fill the entire range of allowed couplings. 
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the branch f—(r) in eq. fl2.4|) . In fact, the AdS vacua associated with the '+' branch all fail 
to satisfy the above inequality. This result can be directly related to the 'ghost' behaviour 
of the gravitons in the latter vacua, which is why we disgarded them in section 0. Eq. ( |3.44| ) 
precisely matches the condition that ensures that the graviton is not a ghost |22j. 

We now turn to the constraints in eqs. ( |3.32| - [3~3"4| ). Using the expression for ti in 
eq. ( |3.31| ), as well as t± = and the definition fl2.5| ) for /oo, we arrive at the following 
constraints: 

(D-3)(D-4)(D 2 -3D + 8) , A . 

Vector: A OB > - <° + >M ~ 3 » , (3.46) 

10 

Sca.ar: A m >_ - <*>^°> ■ (3.47, 

The most stringent constraints here come from the tensor and scalar channels which require 

(3fl- !)(£>- 3) (D-3)(D-4)(D 2 -3D + 8) 

4(^ + 1)2 " Agb " 4(D 2 -5D + W) 2 ^- 48 ^ 

in order that no negative energy fluxes appear in any channel. Similar constraints on the 



GB coupling arise from demanding that the dual CFT is causal [15]. For five-dimensional 
GB gravity, it was found that these two sets of constraints were identical @, Below 
we will show that the flux constraints above again match the causality constraints for an 
arbitrary D. 

4. Gauss-Bonnet black hole perturbations 

The topics of the two subsequent sections will be the study of causality violations (section ||) 
and hydrodynamics (section |6|) in the CFT's dual to GB gravity. In both cases, the analysis 
focuses on examining metric perturbations propagating in the GB black hole background 
( |2.3| ). Hence it is useful to outline the general framework and to derive the equations of 
motion for these perturbations in the present section. To begin, we review the techniques 
introduced in [S3] to study black hole quasinormal modes. The specific background of interest 



is the GB black hole (2.3) introduced above. For the following analysis, it is convenient to 
choose the radial coordinate u = r^/r 2 , in which case the metric becomes 

rl ( , . L 2 



ds 2 = ^-A dt 2 + V dx 2 + -rij^du 2 , (4.1) 

UL \ foe ~t ) ^ 2 f(u) 

with 



l-A/l-4An R (l-n^ 1 



/(«) = " ~ (4-2) 

^GB 

With this radial coordinate, the horizon is now located at u = 1 and the boundary at u = 0. 
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Now on this background (4.1), we wish to study gravitational fluctuations h au . The 
latter are taken to be of the form h uu = h^ u (u)e~ lujt+iqz , where we choose the direction 
z to be the planar coordinate x D ~ 2 . As usual according to the AdS/CFT dictionary, these 
perturbations are dual to stress-energy probes of a finite temperature plasma in the boundary 
CFT. For simplicity, we restrict to the case where the metric perturbations do not couple to 
fluctuations of other background fields. The gravitational fluctuations h^ u can be classified 
according to their transformation properties under the remaining symmetry group SO(D — 3) 
acting in the x % directions — where % = 1, . . . , D — 3, as in the previous section. Hence, htt, 
htz, h zz , h uu , htu an d h zu transform trivially under these rotations and can be considered 
as spin perturbations. 9 Similarly, ht%, h z % and h u % transform as as spin 1 perturbations 
and h%~ transform as spin 2 tensors. Therefore, we have three symmetry channels for gravity 
perturbations, conventionally referred to as the sound, shear and scalar channels, where the 
terminology is adopted from the hydrodynamic description of the dual CFT plasma. The 
equations of motion for fluctuations belonging to different symmetry channels are guaranteed 



to decouple due to the SO(D — 3) symmetry. Following [32|, it is convenient to rescale the 
fluctuations as 

nL 2 /oo - uL \ R - UL \ H-S^H - UL2 — 

(4.3) 

with i,j^t and h = h^. In terms of these, we make a particular choice of perturbations 
for each of the symmetry channels for further study in the following: 



Sound channel (spin 0) 
Shear channel (spin 1) 
Scalar channel (spin 2) 



Htt , Htz , H zz , H uu , Htu , H zu , H, 
Htx , H zx , H ux , 

H xy . (4.4) 



Note that we have renamed x 1 = x and x 2 = y for notational convenience. Using the momen- 
tum and frequency of the fluctuation, we introduce the following dimensionless quantities 

(4-5) 



2vrT' 2vrT 



where T is the Hawking temperature (|2.q ) of the black hole. Furthermore, we can define the 
following "gauge-invariant" variables in each of the channels 

Sound channel: Z sound = q 2 ^ H u + 2n> q H tz + vo 2 H zz + q 2 ( ^ - "j- - % ) H, 

J co \ ioo loo 9 / 

Shear channel: Z shoar = q H tx + to H zx , 

Scalar channel: Z SC ai ar = H xy . (4.6) 



9 Note that the present classification of the probes is identical to that discussed for the states in the scattering 



processes in section 3.3. Unfortunately here, we cannot refer the spin perturbations as scalar perturbations 



because this nomenclature clashes with the hydrodynamic terminology - see below. 
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They are gauge-invariant in the sense that they are invariant under the residual infinitesimal 
diffeomorphisms: h^ u — > — V^ u — V u ^, where £^ = £^(r)e~ iwt+l132 and the covari- 
ant derivatives are taken with respect to the background metric ( |4.1| ), which preserve the 
background metric ( [4.11 ). By first finding the equations of motion obeyed by the rescaled 
fluctuations H^ u , we will be able to derive second-order differential equations for each of the 
Z'b. 

To obtain the equations of motion for H^, we perturb the black hole metric ( |2.3| ) by 
dSp ert , which contains an infinitesimal parameter e. Then, we evaluate the GB Lagrangian as 
a series in e and pick up the coefficient of the second-order term. By varying this resulting 
Lagrangian with respect to each of the fluctuations we obtain the equations of motion of 
interest. We now present each of the three cases above in turn. 

4.1 Scalar channel 

The scalar channel is the simplest case as we only need to consider one perturbation, namely 

r 2 

dsp ert = 2e ^2 ( H xy dx dy) . 
The resulting equation of motion for the perturbation Z scalar = H xy is the following: 

7" _i_ nW 71 _i_ f( 2 ) 7 _n 

scalar ' '-'scalar ^scalar ' '-'scalar "^scalar u 5 

with the following expressions for the coefficients 

r (D _ r 
scalar ufNM ' 

C ( 2 ) _ n 



scalar 



Here and in the following, the expressions denoted by capital calligraphic letters are expres- 
sions that we define in Appendix |C] in order to present our results as succinctly as possible. 
These scalar channel equations were already obtained for GB gravity in general spacetime 
dimensions by [16|. 



4.2 Shear channel 

In the shear channel, we perturb the black hole with 

2 - r2 



dSpert = 2e {Htx dt dx + H zx dz dx + H ux du dx) . 

We choose to work in the gauge where H ux = 0. The resulting equations that we get from 
varying the action with respect to Ht x , H zx and H ux are, respectively 



(D-3) M , {D-l? q-M 
uM tx 16/oo ufM 



H tx o Tr H tx 77T7 71T7 H tx + ™tLzx) - U, 
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H' 7sr H 7~r~r~ H' tx — 0. 



Next, we multiply the first equation by q and the second one by to, add both equations 
and, using the third equation, find a differential equation for the gauge-invariant variable 
^ S hcar = <{H tx + voH zx . The result is 

•^shcar + ^shiar -^shcar + ^shiar "^shoar = 0, (4-7) 

with the following expressions for the coefficients 

AD -2o) 2 / 00 PAA-( J P-3)q 2 / 2 .M 3 

2ufMM{^fM-x y oo M) ' { -> 

4.3 Sound channel 

Lastly, let us consider the following perturbation in the sound channel 
r 2 / f 

ds 2 pcrt = e -^( K j-H tt dt 2 + 2H tz dt dz + (D - 3)H U {dx l f + H zz dz 2 + 

+ 2#£ M du + 2// zu <iu + H uu du 2 ^j . 

Choosing the gauge Ht u = H zu = H uu = 0, the equations of motion that we obtain from vary- 
ing the action with respect to H u , H tz , H = ^ a H aa , H zz , H tu , H zu and H uu , respectively, 
are the following: 

(D - 3) ufM (D - 3) ufM zz 16/oo u/Af ' 

(D-3) A4 , (£»-l) 2 (£>-3) qtx)M 

2 uAA z 16/oo w/AA H " U ' 

H{ t - {D - H" - ^ + ^.H[ t - {D - A)^-^H'+ 
y AA AA u/AA 'u/AA 2 

7> CD-lf q^M (D-l) 2 q 2 X 

fl 2Z " ;^TT7 titz — — J77 ±±tt+ 



ufM 2 zz 8 upM 16/oo u/AT 

~ ^ ^ufW 2 16 upM H "~^ 

H *-( D -*)jf H + u-fN tt ~ { ) ^W~ 2 16 ^PN H ~^ 
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H > i H' zz q (23-1) q/C , (g - 1) g „ 

+ (23 - 3) + (23 - 3)» * z + 2(23 - 3) n/roAA | ^ + 4 n/A/ -§ 

+ 4(2? - 3) u/A ^ H » ~ °' 



q/ g , (D-3)*fM H , (23-1) q/C 

«»/oo " mfooN 4/oo UW J\/i 



V ; z * V * 8(23- 2) V V / /» / 8(23-2) fV 

(23 - l) 2 qW (23-l) 2 qtx>AA 
8(23- 2) /ooV * 4(23-2) /V * 2 

Again, the goal is to use these equations of motion to write a second-order differential 
equation for the gauge-invariant variable 

^sound = -7— H tt + 2wqH tz + t» z 22 22 + — 1 + , 22. (4.10) 

After a series of algebraic manipulations, we arrive at the following equation 

v" _i_ 7' -\- 7 — n fzLi~n 

^souncl "T "-sound ^soimd ~ '-'sound ^sound u ) V 4: -- L - L ,J 

where the expressions for the coefficients cf^d and C^\ ai are extremely long and are given in 
Appendix [p]. 

5. Causality constraints 

Having obtained the equations of motion for the gauge-invariant variables Z in each of the 
channels, we now proceed to use them to study causality violation in GB gravity, following 
[R4L . We will focus on the shear and sound channels because the scalar channel was already 



considered for GB gravity in arbitrary dimensions by [16|. Let us first outline the general 
strategy behind this approach, which will then be used in each of the channels separately. 
First, note that in each of the channels, the differential equation for Z has the form 

Z"{u) + cV>Z'(u) + C (2) Z(n) = 0. (5.1) 

The general idea given in |15[ consists of introducing a new radial coordinate p and rescaling 
the profile Z = ip(u)/B(u) to bring this equation ( |5.l[) into the form of an effective Schrodinger 
equation, 

-d 2 p ip + UTp = m 2 ip. (5.2) 
This form can be achieved with p and B defined by 

dp _ (23 - 1) 
du ~ ~4m 1 /2/( u ) ' 



(5.3) 



D-l 
2 



dlnB 1 (D-l)u- 

= 1 . (5.41 

du 2 Au Auf (1 - 2A GB f ) v ' 
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Further, we introduce h = 1/q and a = to/q. Then upon dividing the above equation (5.2) 
by q 2 and taking the limit q — > oo (or equivalently H — > 0) with a fixed, we arrive at 

-h 2 d 2 t(; + (U + h 2 U 1 + ---)ip = a 2 ?p. (5.5) 

Note that the effective potential in this equation was obtained by Taylor expanding U in the 
limit h — > 0. Note that with the above construction, the potential U° only depends on u, D 
and A GB , whereas U 1 and higher order terms may also depend on a. The key point is that in 
the h — > limit, the dominant contribution to the effective potential comes from U . Hence, 
for our purposes, it is sufficient to study the equation 

-h 2 d 2 rp + U° rp = a 2 ip . (5.6) 

All the preceding is applicable to any of the three channels that we are considering. We 
will now study the behaviour of equation ( |5.6[ ) in each of these channels and give arguments 
as to what conditions have to be satisifed in order to preserve causality in the dual CFT. 

5.1 Shear channel 

The leading potential (as h — > 0) in the shear channel is 

n o [-3 + D - 2{D - 1)A GB - 2(D - 5)A GB / (1 - A GB /)] / 

(^-3)(l-2A GB /) 2 / 0O ' { -> 

while the expression for the subleading potential Ug heat is too long to be presented here. With 
this explicit formula at hand, we can plot the leading potential as a function of u (using ( |2.4| ) 
and ( |2.5D ) in the physical region u € [0, 1]. For any A GB < 1/4, we have that E/° eaI (it = 0) = 1 
and U® heai (u = 1) = 0. For small values of |A GB |, U° is a monotonically decreasing function 
between the boundary and the horizon. However, as we make A GB more negative, the potential 
develops a single maximum between u = and u = 1. This behaviour is illustrated in figure 

I 

The appearance of a new maximum at some < u < 1 implies the existence of quasi- 
normal modes with Re(a 2 ) ~> 1. In turn, this implies that in the limit q — > oo, in which 
we wrote the Schrodinger equation, Re(o>)/q > 1 for these modes, leading to a violation of 
causality in the dual CFT. Hence to avoid causality violation in this channel, we will impose 
a bound on A GB in order to avoid the appearance of a maximum in the potential. 

It is clear from the plot in figure |l| that there is some critical value A GB ar below which U° 
exhibits a maximum. To determine the precise value, we first Taylor expand the potential 
around the boundary u = obtaining 



_ (D - 3) (1 + vra^) + 8A OB ^ 
shcar 2(D - 3) (1 - 4A GB ) 1 ; 

The sign of the coefficient of the n^ -1 ^ 2 term will determine whether or not a maximum 
appears in U° at u > 0. The vanishing of this coefficient then determines the critical value 
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0.2 0.4 0.6 0.8 1.0 



Figure 1: (Colour online) Leading potential U° for the Schrodinger-like equation ( |5.6| ) in the shear 
channel. The blue (solid), green (long dash) and red (small dash) lines correspond to D = 5, 6 and 
7, respectively. For each dimension, the three potentials correspond to Agb = —0.02, —1.5, —3.5, from 
bottom to top. We also show the line U° = 1. The behaviour for higher D is similar. 



A^ = - T • (5-9) 



of the GB coupling: 

(D + i)(Z3-3) 
16 

If the GB coupling is greater than or equal to this value, the potential monotonically decreases 
to at it = 1. Therefore, we obtain the following lower bound on A GB which has to be satisfied 
to preserve causality in the shear channel: 



(£> + l)C?-3) 
16 



Note that for D = 5 and 7, we recover the bounds X^^ 1 = —3/4 and -2 which were originally 
found in |3(| and [[L7j] , respectively. We also observe that eq. ( |5.10D precisely matches the 
bound ( gg§ ) that was derived to avoid the appearance of negative energy fluxes in the vector 



(or spin 1) channel in section 3.S. 



As D becomes very large, we see that A^g" is unbounded from below. This is a clear 
indication that ( |5.10| ) cannot be the correct lower bound for the full GB theory. Indeed, we 
now show that the correct lower bound comes from the analysis in the sound channel. 
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5.2 Sound channel 

The leading potential (as H — )■ 0) in this channel is given by the following expression 



-i 



I. r", , : (D - 2) (1 — 2A GB /) 2 foo(D — 3 — 2(D - 1)A GB - 2(D - 5)A GB / (1 - A GB /)) 
(D - 2){D - 3) - 6(D - l)(D - 2)A GB + 12(D - 1) 2 A GB + 

- 2A GB / (1 - A GB /) ((D - 2)(D - 9) + 12(D - 1)A GB + 

- 2(D -3)(D- 5)A GB / (1 - A GB /))] / . (5.11) 

The behaviour of this function is the same as that described in the shear channel, as can be 
seen in figure [2|. 

if 




Figure 2: (Colour online) Leading potential U° for the Schrodinger-like equation in the sound chan- 
nel. The blue (solid), green (long dash) and red (small dash) lines correspond to D = 5, 6 and 7, 
respectively. For each dimension, the three potentials correspond to A GB ~ —0.02,-1.5,-3.5, from 
bottom to top. The behaviour for higher D is similar. 

Following the same steps as in the previous case, we obtain the following critical value 
for the GB coupling 

, sound _ (£>-3)(3£>-l) 



A 



CD 



4(D + 1) 2 ' 

Hence, by the same arguments used before to avoid causality violation in the theory, we get 
the following lower bound 



Av. \ sound 
GB ^ ^ GB 



(D-3)(3D- 1) 
4(D + 1) 2 ' 



(5.12) 
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Again, for D = 5 and 7, we recover A^g 1 " 1 = —7/36 and -5/16, as was originally obtained in 
[31, 30 1 and jl7|], respectively. Again comparing to the results in section 3.3, we also find that 
eq. (5.12) precisely matches the bound (3.47) necessary to avoid the appearance of negative 
energy fluxes in (what was denoted there as) the scalar or spin channel. 

Eq. ( |5.12| ) provides a more stringent lower bound than that in the shear channel ( |5,10|) . 



In particular, when D becomes very large, A^g 1 " 1 is bounded from below by —3/4. 
5.3 Scalar channel 

For completeness, we present the analysis for the scalar channel, which was already considered 



in 1 16]. The expression for the leading potential is in this case 



U~ 



(D - 4) (1 - 2A GB jf /«, (D - 3 - 2(23 - 1)A GB - 2(D - 5)A GB /(1 - 
(D - 3)(D - 4) - 2(D - 1)(D - 6)A GB - A(D - 1) 2 A GB + 
- 2A GB / (1 - A GB /) [3 (D(D - 7) + 14) - 4(23 - 1)(223 - 7)A GB + 



Agb/)) 



2(23-5)(23-7)A GB /(l-A GB /)] 



(5.13) 



The behaviour of this function is shown in figure |3], where we have made use of the bounds 
found in the shear and sound channel to choose the minimum values that A GB can take. 

if 




Figure 3: (Colour online) Leading potential U° for the Schrodinger-like equation in the scalar chan- 
nel. The blue (solid), green (long dash) and red (small dash) lines correspond to D = 5, 6 and 7, 
respectively. For value of D, the three potentials correspond to A GB = —0.05,0.15,0.21, from bottom 
to top. 

We note that, as opposed to the two previous cases, here the single maximum in the 
potential appears when we increase the value of A GB . The critical value of the coupling in 
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this channel is given by 

(D - 3)(D - A)(D 2 - 3D + . 



X s 



4(D 2 - 5L> + 10) 2 

Hence, in this case this value gives an upper bound on the GB coupling 

GB " 4(jD2 _ bD+W) 2 • I 5 " 14 ) 



For D = 5 and 7, we recover the results Aq^ t = 9/100 and 3/16, first obtained in [15[ and 



[17 1, respectively. Eq. ([5.12D also precisely matches the upper bound (|gD necessary to avoid 



the appearance of negative energy fluxes in the tensor or (spin 2) channel. Further, we see 
that as D becomes very large, Aqj 1 " approaches the critical value A G b = 1/4 from below. 

5.4 Plasma instabilities 

To complete this section, let us make some comments on other possible instabilities in the 
theory. So far, our analysis of the effective Schrodinger equation (|5.6|) concerned the possible 
appearance of superluminal signals in the dual CFT. However, it was noted in [jl4| that in 
D = 5 GB gravity, a new instability arises in the scalar channel at A GB = — 1/8. If one 
decreases the coupling below this value, the potential i/£ alar develops a negative minimum 
located right in front of the horizon u = 1. Following [pSfl , it was argued that for large q, this 
well supports unstable quasinormal modes. This effect was extensively studied for charged 
black holes in GB gravity in [16, [34J. In particular, the results of [fill] indicate that a similar 
well develops in the scalar channel of GB gravity for general D, e.g., see figure ^ To obtain 
the value at which this new instability occurs, expand Q5,13| ) to around u = 1 

uQ _ (1 + VI - 4A GB ) {(D - 3)(£> - 4) - 2(D - 1){D - 6)A GB - 4(D - 1) 2 A GB ) 



scalar 



D-l 

u~ -11 + 



2{D-A) (D - 3 - 2{D - 1)A GB ) 

The sign of the coefficient of (u^ -1 )/ 2 — l) will determine whether the potential develops a 
negative well in front of the horizon. One finds that the negative well develops for [16 

(Z3-6) + V5D(Z3-8) + 84 
Agb < . (5.15) 

Comparing this expression with the lower bound on A GB obtained in the sound channel ( |5.12| ), 
we find that this instability arises within the well-behaved regime only for D = 5 and 6. That 
is, for 

7 1 

D = 5: < A GB < — , 

36 ~ 8 

D = 6: - iM <A OB <-J- (5.16) 
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the GB theory is dual to a well-behaved CFT, but the plasma still seems to exhibit instabilities 
corresponding to unstable quasinormal modes in the black hole background. Given that these 
modes have large momentum, it seems that the instability causes the homogeneous CFT 
plasma to "clump" into an inhomogeneous configuration. However, for D > 7 (or d > 6), the 
values given by ( |5.15| ) are ruled out by the lower bound ( |5.12|) and so no such scalar channel 
instabilities arise in the dual CFT's (which are not already pathological for other reasons). 




Figure 4: (Colour online) Instabilities in C/ S ° calar . The blue (solid) lines correspond to D — 5 for 
A G b = —0.15,-0.19, whereas the red (dashed) lines correspond to D = 6 for A GB = —0.245,-0.26 
(the well is barely noticeable in this case). 

Of course, one should also examine the shear and sound channels for similar instabilities. 
In [31], it was observed for D = 5 that the shear channel does not produce any new instabil- 



ities. Similarly, examining the potential (5.7) the same result applies for arbitrary values of 
D. It was also observed in j 31 ] that in D = 5 the sound channel can develop instabilities for 
A GB > 1/8. Hence, this effect only occurs outside of the physical regime, i.e., A GB < 9/100. 



For D > 6 we examine the potential (5.11) for the possible appearence of a negative well. We 
find U® ound > everywhere in the range 1 > u > for D > 6. Hence if we restrict the GB 
coupling, as in eq. ( |3.48 ), to produce a physically well-behaved CFT, the shear and sound 
channels do not produce any new instabilities for a homogeneous plasma. 



6. Second-order hydrodynamics 

Hydrodynamics of relativistic fluids is discussed in |35j. Here we summarize a few relevant 
results for our computations. Let us consider a (D — l)-dimensional conformal fluid, as 
is appropriate for the CFT dual of D- dimensional GB gravity. The dynamics of the long 
wavelength fluctuations can be organized in terms of a derivative expansion. In particular, 
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in the absence of any conserved charges, the dynamics of the hydrodynamic modes is simply 
governed by the conservation of the stress-energy tensor: \/ a T ab = 0. The latter includes both 
the equilibrium part, involving the local energy density e and pressure P, and a dissipative 
part U ab 

T ab = e u a u b + PA ab + U ab , (6.1) 

where A ab = g ab + u a u b and u a is the local four- velocity of the fluid, with u a u a = — 1. For a 
conformal fluid, we also have the restriction that the trace of the stress tensor must vanish, 
i.e., T a a = 0. In the equilibrium contribution above, this requires that P = e/(D — 2). 

The dissipative term can be written as an infinite series expansion in velocity gradients 
and curvatures (for a fluid in a curved background), where the coefficients in this expansion 
are the transport coefficients of the fluid. At first order for a conformal fluid, one has 



where 



nf = — Tj <j ab . 



a ab = 2V (a u b) = A ac A bd ^ ^ + y ^ _ ^ ^ 



Implicitly here, we have used the condition that T a a = to eliminate any bulk viscosity 
contribution. Hence the only nonvanishing transport coefficient appearing at this order is the 
shear viscosity rj. Of course, if we truncate the hydrodynamic equations to include only this 



"Newtonian" first-order term, there are superluminal modes propagating in the fluid [36]. 



One can try to avoid this problem by going to next order in the expansion [37]. In general 
for a conformal fluid, five new transport coefficients will appear in the second-order term 
nf. However, there is only one contribution which is independent of both the vorticity and 
background curvature, as well as linear, 



nf = r?r n 



{ u ■ V<x ob > + ~^f ah (V • it) 



The additional transport coefficient above is the relaxation time rn and this linear term is 
sufficient to examine the question of causality within the hydrodynamic regime . 

In J|], the authors studied linearized fluctuations of the second-order hydrodynamics of 
conformal fluids and found the following results: 

• Shear channel 

The dispersion relation in the shear channel was found to be 

9 „ IVO nil 

-tx> 2 rnT- — +q 2 ^ = 0, 6.2 
2ir s 



where tx> and q are defined in eq. (|4.5| ). We are interested in the hydrodynamic limit in which 
vo, q — > 0, with t»/q kept fixed. In this limit, the dispersion relation ( |6.2D yields the following 
Taylor series solution 

2 

vo = -2vri 5 q 2 _ 8vr 3 i r n T % q 4 + Q(q 6 ). (6.3) 
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Note that we have discarded a solution where to remains finite as q — > 0, hence lying beyond 
the hydrodynamic regime — see j| . The wave- front speed with which disturbances propagate 
out from a discontinuity in any initial data is governed by [^] 

C S| h»*4.jX. (6.4) 

Thus, requiring causality to be preserved in this channel of a conformal fluid imposes the 
following bound [31] 

r u T > 11 . (6.5) 
s 

• Sound channel 

The dispersion relation in the sound channel was found to be 

* m i™ 2 9 9 m ?2(Z) — 3) n iq 2 c 2 
27T {V — 2) s 2ir 

where c s is the speed of sound, which in any conformal fluid is a constant 

At small momentum q, the Taylor series solution to the above equation corresponding to the 
sound wave is 

to = c s q - 2vri TT q 2 + — rT f c 2 r n r - — J q 3 + 0(q 4 ), (6.8) 

where 

IT-§Z|n. (6,) 

Note that we have only written one of the three solutions to ( |6.6| ). The second solution 
corresponds to waves propagating in the opposite direction, i.e., replace q with — q in eq. ( |6.8| ). 
The third solution again lies beyond the regime of validity of hydrodynamics. Given eq. 
we have 

f * Re(vo) r. 2(L>-3)r? 1 

vf:Z = t ™ — H = \ c + m o\ ^r- 6 - 10 
|q|^oo q y [D -2) s r n T 

Hence, causality in this channel imposes the following bound 

r n T>2-. (6.11) 

s 

It is quite interesting that despite the explicit appearance of D in eq. floTToD , the bound ( jsTTH ) 
is independent of the dimension and takes the same simple form as originally found for D = 5 
in §![]. 

We see that in both channels considered above, the wave-front speed of the linearized 
fluctuations diverges as the relaxation time rn goes to zero. Notice also that the front velocity 
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in the sound channel is always larger than that in the shear channel. Therefore, the bound 
( |6.11 ) provides a stronger constraint on the transport coefficients. As emphasized in [31|, the 



constraints ([O^) and ( |6.11| ) arise when considering linearized modes outside the regime of 
validity of hydrodynamics. Hence, these constraints must not be regarded as fundamental, as 
they merely indicate where a certain approximate framework describing the conformal fluid 
becomes problematic. 

6.1 Holographic Gauss-Bonnet hydrodynamics 

Next we would like to apply the now standard techniques [4C , f|] of the AdS /CFT framework to 
calculate the transport coefficients rj and Tn for the CFT's dual to GB gravity in an arbitrary 
spacetime dimension. The shear viscosity for GB gravity in D = 5 was originally calculated in 



[14 1, however, these calculations are easily extended to GB gravity in any dimension [14, 16]. 



Their result is 



rj 1 

s 4-7T 



1 2{D ~ l) \ 
1 " {D _ 3) 



.12) 



where the entropy density s is given in eq. ( p.8| ). We reproduce this result with our analysis 
of the sound waves below. 

Hence it remains for us to calculate the relaxation time tr. For this purpose, we must 
solve the equation of motion ( |4.11| ), which describes the propagation of sound waves in the 
CFT plasma dual to GB gravity. This allows us to obtain expressions that can be compared 
to those presented in the previous subsection. First, we make the following ansatz for the 
solution near the horizon 

^sound = (1 — It 2 ) , 

where there could be an overall coefficient independent of u, which is not relevant for our 
discussion. If we replace this ansatz into ( [4.11 ), we readily obtain the following two solutions 



with the negative sign satisfying the correct infalling boundary condition at the horizon. Let 
us now focus on the hydrodynamic limit, where to, q — > with ™ kept fixed. We will write 
our solution as a series in q and then solve eq. ( [4.11 ) perturbatively. The horizon boundary 



condition implies that we can write our solution as 

Z sound = (l-u^) - ^ (z (n;t»,q) +iqzi(n;t»,q) +q 2 z 2 (n;t»,q) + C(q 3 )) , (6.13) 

where at the horizon lim M _»i Zi(u;to, q) = 6® with i = 0, 1,2. Further each of the Zj(u;t», q) 
should only depend on the ratio t»/q which implies a scaling invariance 

Zi(u; fia, ficj) = Zi(u;xx>,cfi (6.14) 

for any [i. In this expansion, we have only written terms to the order necessary to identify the 
sound wave dispersion relation fl6.8|) , which is obtained by imposing the Dirichlet boundary 



- 25 - 



condition on Z sonnd at the asymptotic boundary 

lim Z sound = . (6.15) 

■u-i-0 

By comparing the terms of the dispersion relation at each order in q, we will be able to identify 
the expressions for the conformal sound speed c s (the leading order), the shear viscosity to 
entropy ratio rj/s (the first order) and the relaxation time rn (the second order) . The analysis 
is greatly simplified by adopting a new radial coordinate 



l-4A GB (l-uV). (6.16) 
With this radial coordinate, the horizon and asymptotic boundary conditions become 

lim Zi(x; vo, q) = 6^ (horizon), 

X-— s-l 

lim ^sound = (boundary). (6-17) 

6.1.1 Speed of sound 



Solving eq. ( 4.11 ) perturbatively, the leading order term gives a second-order differential 
equation for zq{x). We will not write the explicit equation here, but merely state its solution, 

= q 2 [{D - 5)x 2 + 4x - (D - 1) (1 - 4A GB )] - 8 (D - 2) x 
Z ° 4[( J D-l)q2-2( J D-2)/ 00 tx,2]A GB x ' ( ' ' 

By imposing the Dirichlet boundary condition ( |6.17 ) and solving for vo, we find the expected 
dispersion relation 

1 

ro = , q . 6.19 



Of course, the factor — 2 is precisely the speed of sound (6/7) for a conformal fluid in 

D — 1 dimensions. 



6.1.2 Shear viscosity-to-entropy ratio 

To first order in q we obtain a second-order differential equation for z\ (x) , which involves the 
previously found zq{x). This equation is even more involved than that for zq{x). To proceed 
we make the following ansatz for the solution z\ = (t»/q) zq{x) F(x). Again, it is not very 
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illuminating to write the resulting equation and so we simply give its solution 

F{X) = q 2 [(D - 5)z 2 + 4x - (D - 1) (1 - 4A GB )] -8(D-2) /octx> 2 A GB x * 

[q 2 [(D - 5)x 2 + 4x - (D - 1) (1 - 4A GB )] - 8 (D - 2) joon> 2 A GB x] x 



(L>-1)(1-4A GB ) 



4(L> - 5) 



arctan 



V(£> 5) (1-4A 



( J D-l)(l-4A GB ) + (D-5)a; 



(1-x) + 



2 V 2 



(£> - 1)(1 - 4A GB ) 5)x 2 + (£> - 1)(1 - 4A GB ; 



4(Z>-5) 



2(D - 3) - 4(D - 1)Aqb 



+ 



4q 2 x(l - x) ((£> - 3) - 2(D - 1)A GB ) 



.20) 



Note that while at first sight it might seem that this result is not well-behaved at D = 5, 
taking the limit D — > 5 carefully, recovers the results obtained in |3lJ — in taking this limit, 
one should recall that A GB < 1/4. Given this general solution ( |6.20| ), we can impose the 
Dirichlet boundary condition at infinity and solve for t», obtaining the following dispersion 
relation up to second-order in q 



■ (P-3) 
'2(D-2) 



1 



2(^-1) 
(D-3) 



A 



GB 



q 2 . 



By comparing to eqs. (6.5) and (|6.9| ), we see that the factor in brackets in the second-order 
coefficient is essentially the ratio of the shear viscosity to entropy density with 



V 
s 



1 

47T 



(D-3) 



.21) 



Of course, this precisely reproduces the result found previously in [14, 16 1. 



6.1.3 Relaxation time and causality violation 

To second order in q we obtain a differential equation for 22(2;), involving zq (x) and zi(x). 
Unfortunately, we were only able to solve this equation numerically. 10 Naturally, this means 
that we are not able to obtain general expressions for arbitrary values of D. Here we present 
numerical results for D = 5, 6, 7, 8 but these are indicative of the behaviour for higher dimen- 
sions as well. Note that the constraint ( |6.11| ) to avoid causality violations requires 

t u T-2- > 0. 
s 

We show the results of our numerical computations for this quantity in figure ||. We see that 
for D > 6, causality does not impose an upper bound on the GB coupling, however, we still 



10 The details of numerical calculations will not be presented here, but are available from the authors upon 
request. 
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require A GB < 1/4 as was assumed from the very beginning of the discussion. From these 
numerical results, we find the following bounds for having causal second-order hydrodynamics, 

D = 5 : -0.711 < A GB < 0.113, 

D = 6 : -1.7 < A GB < 0.171, 

D = 7 : -3.1 < A GB < 0.25, 

D = 8 : -5.0 < A GB < 0.25. (6.22) 

Note that we do not quote the bounds in higher dimensions with the same accuracy because 
the numerical analysis became more unstable for large values of D. In figure |(| we also present 
the behaviour of the front velocity ( |6.10D in the sound channel for these various cases. 

t„T-2± 

s 

0.35 1 




Figure 5: (Colour online) Causality in the second-order GB hydrodynamics is violated once the value 
of TuT — 2^ becomes negative. From bottom to top, the blue, green, red and purple curves correspond 
to D = 5, 6, 7 and 8, respectively. 



7. Discussion 

In this present paper, we have made a broad study of holographic GB gravity in arbitrary di- 
mensions. The first step was to construct the AdS/CFT dictionary which applies for arbitrary 
D in section |3| There are two dimensionless parameters which characterize GB gravity, i.e., 
the ratio of the AdS curvature scale to Planck scale, L/£ P , and the GB coupling, A GB . Since 
we are only dealing with the gravitational sector of the AdS theory, these parameters are 
naturally related to CFT parameters which determine the behaviour of the n-point functions 



of the stress-energy tensor. So for example, using eqs. ( 3.15 ) and ( 3.31 ), we can translate the 
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-5 -4 -3 -2 -1 CB 0.00 0.05 0.10 0.15 0.20 0.25 CB 

Figure 6: (Colour online) Front velocity in the sound channel. From top to bottom, the blue, green, 
red and purple curves correspond to D = 5, 6, 7 and 8, respectively. The dashed vertical lines indicate 
where this velocity reaches one on the various curves. 



gravitational parameters to expressions in terms of the central charge Ct and the scattering 
parameter ti. Alternatively, they could be expressed in terms of the three-point couplings 
A, B, C using eqs. ( 



|) and fl3,9|) (keeping in mind the constraint ( 3.10 )). Given the results of 
section [3|, one can express any physical quantities, e.g., tj/s, entirely in terms of parameters 
appearing in the CFT. This is, of course, an essential step in the program of using holographic 
models to characterize the properties of the quark-gluon plasma [41]. 

Using these results, an interesting observation can be made with regards to the relation 
between the central charge Ct appearing in the two-point function and the "central charge" 
characterizing the entropy density. That is, one can use the entropy density of a CFT in 
d dimensions to provide an alternative definition of a central charge: Cs oc s/T d ~ 1 . In 
two dimensions, Ct and Cs are related by a simple numerical factor. However, there is no 
evidence of any simple relation for strongly coupled CFT's in higher dimensions, e.g., see 



[42, 43 1 . This alternative central charge was recently considered for CFT's dual to Einstein 



gravity in [44] and here we extend the discussion to consider GB gravity in an arbitrary 
number of dimensions. First we adopt the normalization of [44] in defining Cs as 



C S 



d + i ( d rr((d + i)/2) s 



d-l\ 2^ 3 / 2 



7T 



jid-l 



(7.1) 



Then examining the entropy density (2^) for GB gravity along with the holographic expres- 
sion for Ct in eq. ( |3.15 ), we see that these two central charges have a rather complicated 
relation: 



Cs = Ct 



{l + 2 ld t 2 ) a 

\d-l 



(7.2) 



;i+7^ 2 ) c 

where "fa = (d — 2){d — 3)/(4d(d — 1)). The normalization was chosen above so that with 
*2 = {i.e., for the CFT dual of Einstein gravity) we have Cs = Ct- In general, one sees that 
Cs is proportional to Ct, however, it also has a highly nontrivial dependence on t-z- Hence we 
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are lead to conclude that these two central charges are independent parameters in the family 
of CFT's dual to GB gravity. 

We may go further by recalling that t 2 is constrained by eqs. ( 3,32 - 3.34 ). With t& = 0, 
combining the tensor ( 3.32| ) and scalar ( |3.34 ) constraints yields 

d- 1 



d-3 



<t 2 <d-l 



(7.3) 



We may note that in this range ( |7.3| ), the ratio Cs/Ct is a monotonically increasing function 
of t% for any dimension. Hence these constraints in turn lead to 



3d + 2 (2{d + 2 



4d 



C s d 2 
< — < — 

~ C T ~ 4:d 



d + 6 (2{d 2 - 3d + 6) 
d 2 - d + 6 



3d + 2 

Explicitly evaluating these results ( |7.4[) for various dimensions, we find 

d = 4: 0.4723 < Cs/Ct < 1.7147 , 



(7.4) 



d = 5 
d = 6 
d = 7 



0.3220 < Cs/Ct < 3.6714, 
0.2185 < Cs/Ct < 8.4280, 
0.1477 < Cs/Ct < 19.6317. 



(7.5) 



Here we see that the maximum possible value of Cs/Ct grows monotonically as d increases 
while the minimum ratio is monotonically decreasing and asymptotically approaches zero. 

Following [44 1, we can compare these results for the strongly coupled holographic field 
theories to those for free fields. First, substituting the free field results (|3.39| ) for A, B and C 
into eq. ( |3.5| ) yields 



Ct 



1 



d 



d 



d l 



d 



- n s + - n f + — n t 



(7.6) 



where n s , fit and fit denote the number of (massless) degrees of freedom contributed by scalar, 
fermion and tensor fields, respectively, as described after eq. ( |3.39| ). It is straightforward to 
calculate the entropy density for free massless fields, e.g., see |f43[| . Then with the definition 
(PD, we find 

d 



C s = 2d(d+1) (JL\ T(d-l)((d) 



n s + 1 



tl-d 



hf + ht 



(7.7) 



Combining these results, we arrive at 



Cs 
Cr 



free 



7T 



(d 2 - 1) 



where 



K{d) 

K{d) 



1 + 



\2ir 
d-3 



r((d-l)/2) n s + (l 



2 l ~ d ) hj + ht 



T{d/2) 



2{d-l) 



t 2 + 



1 



n s + ^hf + 



d -^h t 



4(1 -2~ d ) d 2 -l 



16 



d y r((d- 1)/2) 



2tt 



r(d/2) 



C(d). 



(7.8) 
(7.9) 
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For example, with d = 4, this free field result becomes [44] 



Ct 



8 n s + lhf + 2n v 4/ 1 2 

= 2— 1 = - 1 + -t 2 H 1 4 

free,d=A 3n s + ±hf + l2n v 3\ 6 15 



(7.10) 



Note that in this case, the pre-factor in eq. ( |7.8| ) yields 1/K(d = 4) = 3/4, which corresponds 
to the ratio of the entropy density of N = 4 super- Yang-Mills at strong and at weak coupling 

If we focus on ^4 = 0, as for the dual of GB gravity, eq. (|7.8|) reduces to the rather simple 



expression 



Ct 



K{d) 



free 



1 + 



d- 3 



*2 



(7.11) 



which we might compare to the analogous result ( |7,2|) at strong coupling. The constraints 



( 3.32 - 3.34 ) , derived from considerations of the energy flux, make no reference to the strength 
of the coupling and so eq. (|7.3j) applies equally well for the free field theories. Hence this 
constraint on t 2 in turn leads to 



< 



d- 1 



K{d) 



(7.12) 



free 



Explicitly evaluating these results ( |7.12| ) for various dimensions, we find 



d = 


4 : 


0.6667 < C s /C T \ free 


< 


2. 


d = 


5 : 


1.0884 < C s /C T \ free 


< 


4.3536, 


d = 


6 : 


2.2781 < C s /C T \ free 


< 


11.3906 


d = 


7: 


5.7890 < C s /C T \ free 


< 


34.7338 



(7.13) 

which are readily compared with the corresponding results at strong coupling in eq. ( |7.5| ). 
Here we see that both the maximum and minimum possible value of Cs/Ct |^ ree are growing 
monotonically as d increases. Further the ranges of Cs/Ct for the free fields and the strongly 
coupled holographic CFT partially overlap, however, generally for a given value of t 2 , the free 
field result is the larger of the two ratios. 

We must emphasize, however, that the AdS/CFT dictionary of section [| is by no means 
complete. Eqs. ( |3.15| ) and (3.31) establish a connection between the two gravitational param- 
eters and two CFT couplings characterizing the two- and three-point functions. However, the 
GB action ( |2.1| ) fixes the form not just of these CFT correlators but also an infinite series of 
n-point functions for the stress tensor. Implicitly then, the parameters controlling the higher 
n-point function are not independent couplings within the family of CFT's dual to GB gravity 
or any simple gravitational action with a finite number of couplings. 

Among the interesting results found here were the bounds on the GB coupling in eq. ( |3.48| ) . 
In section ||, these constraints were derived by demanding that the energy flux remained pos- 
itive everywhere for certain scattering experiments in the dual CFT. However, in section 0, 
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precisely the same constraints were reproduced by demanding that causality was respected 
in the dual CFT. We should also note that while it did not set a fundamental constraint, 
causality violations and negative energy fluxes also appear in the shear channel at precisely 
the same critical value of A GB given in eq. fl5.9| ) . The exact match of the constraints derived in 
these two approaches was already observed for GB gravity in D = 5 and 7 in [|30], and |l7j , 
respectively. It was argued in ]3(| that this precise matching is a consequence of the special 
two- derivative nature of the GB gravity. While, as explicitly seen in appendix Hi the energy 
flux calculations only depend on the three-point function of the stress tensor, in general, one 
should expect that the causal propagation of signals in, e.g., a CFT plasma is also effected 
by higher n-point functions. Hence in general, demanding causal propagation in the CFT 
provides independent constraints on the parameters of the dual gravitational parameters. 
Explicitly evaluating the constraints in eq. ( |3.4£ ) for various dimensions, we find 



D = 5 
D = 6 
D = 7 
D = 8 



-0.1944 < A GB < 0.0900, 
-0.2602 < A GB < 0.1523, 
-0.3125 < A GB < 0.1875, 
-0.3549 < A GB < 0.2076. 



(7.14) 



One can compare these results to the constraints (6.22) arising from demanding sublumi- 
nal propagation of sound waves in second-order hydrodynamics. As expected, the above 
constraints are always more stringent than those coming from the analysis of second-order 
hydrodynamics, as was observed in [|3l| for D = 5. However, as emphasized there, these two 
sets of constraints stand on a completely different footing. If the GB theory lies outside of the 
range specified by the causality or energy flux constraints ( |3.48| ) , then the gravitational theory 



and its dual CFT are fundamentally pathological. In contrast, the constraints (6.22) arising 
from second-order hydrodynamics only indicate when a certain approximate description of 
the CFT plasma becomes problematic. 



One can combine the fundamental constraints (3.48) with the hydrodynamic analysis to 
examine the behaviour of the limiting wave-front velocities in the shear ( |6.4| ) and sound ( |6.10| ) 
channels. In figure 0, we exhibit the wave- front speeds in the CFT plasma as a function of 
d = D — 1 for GB gravity with the lower and upper values of A G b allowed by the bounds 
( |3.4§ ), as well as for Einstein gravity (i.e., A GB = 0) — compare to figure 1 of f35]| . Figure 
H also shows the detailed behaviour of the wave-front velocity in the sound channel as a 
function of A GB for each dimension. The general trend is that the wave-front velocities are 
monotonically decreasing as d grows. Given the analytic results for Einstein gravity [p5fl , one 
sees that in this case for large values of d: v^°^ v ~ \/l/d and v^°^ d ~ \/3/d. As seen in figure 
[<], for d > 8, the GB gravity with the largest (smallest) allowed value of A GB has the smallest 
(largest) wave-front velocities. 

By combining eqs. (3.48) and ( |6.21| ), one obtains bounds on the ratio of the the shear 
viscosity to entropy density in the CFT's dual to GB gravity. In particular, there is a lower 
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Figure 7: (Colour online) Limiting wave-front speeds for the shear (bottom) and sound (top) channel 
fluctuations in conformal second-order hydrodynamics. The black bullets correspond to the values for 
Einstein gravity {i.e., Agb = 0). The red squares and blue triangles represent the wave-front velocity 
for the lower and upper values of A GB allowed by the fundamental bounds (3.- 



bound which depends on the dimension of the bulk spacetime as: 



1 > — 
s ~ 4tt 



1 



(D -l)(D- 4) (D 2 -3D + 8) 
2 (D 2 -5D + 10) 2 



(7.15) 



For D = 5, we recover the result rj/s > 16/25 x (1/4tt) = (0.640)/(4-7r) originally found in [15]. 
Recently, it was observed [17| that this function ( |7.15| ) has a minimum in the vicinity of D = 9, 
as illustrated in figure ||. More precisely, one finds a minimum value of r\j s ~ (0.414)/(47r) for 
D ~ 9.207. Further for large D, this lower bound rises again and asymptotically approaches 
rj/s = 1/ (87r) . It seems that these results may provide interesting empirical data in considering 
the question of a fundamental lower bound for rj/s, as was conjectured by j/J. However, it is 
not clear to us what conclusions one may draw with respect to this question at this point. 

Within this class of CFT's dual to GB gravity, one also finds an upper bound for the 
ratio of the shear viscosity to entropy density 



< — 

S 47T 



1 + 



(D — 1)(3.D — 1) 
2(D + 1) 2 



(7.16) 



As shown in figure |9|, this upper bound is a monotonically increasing function which ap- 
proaches r//s = 5/(87r) asymptotically for large D. Of course, there is no fundamental in- 
terpretation of this upper bound since r\js can become arbitrarily large for weakly coupled 
theories. Again, this bound only applies for the strongly coupled CFT's with a holographic 
description in terms of GB gravity. 
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Figure 8: (Colour online) Wave-front speed in the sound channel in conformal second-order hydro- 
dynamics in d — 4, 5, 6, 7, 8, 9 (from top to bottom). The dashed vertical lines indicate the upper and 
lower values for A G b allowed by the bounds ( 3.48 ). We see that for d = 8, 9, the speed for the upper 
value of Aqb is lower than the Einstein result Agb = 0. This is in agreement with the behaviour shown 
in figure J?]. 





Figure 9: (Colour online) On the left, we plot the lower bound ( 7.15 ) on the ratio rj/s as determined 
by the upper bound on A G b imposed by requiring a causal CFT for various dimensions d = D — 1. 



On the right, we plot the analogous upper bound ( |7.16|) on the ratio r]/s within the class of CFT's 
described by GB gravity. Note that for d — 4 and 5, the bullets indicate the upper bound taking into 
account the plasma instabilities discussed in section |5.4|. 



As also shown in figure ||, the previous upper bound is probably slightly lower in d = 4 
and 5 when one takes into account the plasma instabilities discussed in section 5.4 . These 



instabilities do not indicate a fundamental pathology for the GB theories in the ranges given 
in eq. ( |5.16j ). Rather in these cases, one simply concludes that a homogenous plasma is an 
unstable configuration in the dual CFT. Since the unstable modes have large momentum, it 
seems that the plasma wants to "clump" into an inhomogeneous configuration. The immedi- 
ate implication for our analysis is that the hydrodynamic calculations, and in particular, the 
computation of rj/s, is unreliable in this regime given in eq. ( |5. 16 ). From a gravitational per- 
spective, the plasma instabilities correspond to unstable quasinormal modes of the black hole 
background. It would be interesting to carry out a more extensive analysis of the quasinormal 
mode spectrum in GB gravity. This would likely provide greater insight into the nature of 
these instabilities. 
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A. Energy fluxes in terms of the three-point couplings 



In this appendix we present an explicit computation of the energy flux (|3.6| ) for a CFT in 
an arbitrary number of dimensions. In particular, t<i and £4 are determined in terms of the 



coefficients which appear in the three-point function of the stress-energy tensor [27, 26 1. Our 
calculations closely follow those given in [17]. 



Let us consider the expression in eq. (3.16). Without loss of generality we can position 
the detector on the (positive) x 1 axis, i.e., n l = 8 l \ as in section |3.3| . We then introduce the 
light-cone coordinates = iix 1 and express the energy flux measured at large distances as 



lim 

x+— ¥00 



d-2 



dx T (x 



,x 



(A.l) 



In order to determine ti and £4, it suffices to compute the energy one-point function ( |3.16 ) with 
specific polarizations which yield two independent linear combinations of these coefficients. 



Recalling the discussion at the beginning of section 3.3, we will consider polarizations in the 
tensor and vector channels. For the tensor channel, we choose a polarization where the only 
nonvanishing components are ejj = ejj where 1 and j are fixed indices with i, j > 1 and 1^3. 
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Similarly for the vector channel, our polarization will only have nonvanishing components 
e n = e n where i is again a fixed index with i > 1. The corresponding linear combinations of 
^2,^4 which will appear in eq. ( |3.6[ ) are those appearing in the constraints ( p.32|) and ( |3.33| ). 
That is, 



Tensor : {S(n)) 
Vector : {S(n)) 



E 



E 



1 



1 



d-1 



t-2 



d 2 - 1 



a, 



d-3 



d 2 - 1 



u 



(A.2) 
(A.3) 



We will present the details of the calculations for the tensor channel below. The vector 
channel calculations are a straightforward extension of these. 

So choosing the tensor polarization above, the numerator of eq. ( |3.16| ) is proportional to 
the following three-point function 



/: 



,(*> = / 



d d x e iEt lim 



1 Xi 



d-2 



dx 1 { T V3 (x) T — (xi ) T Vj (0) ) 



(A.4) 



(where again i and j are fixed indices and so no sum is intended above). Similarly, the 
normalisation in the denominator is provided by the two-point function 



h{E) 



d d x e iEt (T i3 (x)TdO)}. 



(A.5) 



In what follows we find it convenient to assume that the total number of space-time 
dimensions is even. In particular, this assumption allows to use the residue theorem in 
the evaluation of certain integrals. Since our final result turns out to be insensitive to the 
parity of the spacetime dimension, we analytically continue it (and implicitly the intermediate 



integrals) to any d. Let us proceed with the evaluation of the numerator ( A.4 ) and the 
denominator ( |A.5| ) separately. 

We start with the denominator QA.5| ) where we can simply use the general result for two- 
point function ( |3.1| ). Inserting the indices i] and explicitly evaluating the general expression 
yields 



(T»(x)T„(0)) 



Ct 



2(x 2 



1 



((x ? ) 2 + (xi) 2 ) + ^(x J ) 2 (xi) 2 



The central charge Ct is given in eq. (|3.5|), which we present here again 



Ct 



d 
7T2 



(d - l)(d+ 2)A -2B- 4(d + 1)C 



r(d/2) 



d(d + 2) 



(A.6) 



(A.7) 



The constants A, B,C are the parameters which control the three-point function of the stress 



tensor [27, 26 1. To evaluate the integral in eq. (A.5), we must first provide an ie prescription for 
the operators. Recall that the operators Tjj in this expression are not time-ordered which then 
requires that t—tt — ie and so the light-cone coordinates are both replaced by x 



± 



ie. 
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The various integrations are then best performed by first integrating over spatial directions 
which are perpendicular to x , x % and x 3 using (d— 4)-dimensional spherical polar coordinates, 
subsequently integrating over (x*,x^)-plane using polar coordinates and finally the integrals 
are performed using the residue theorem for the poles in these expressions. The latter integrals 
are 



dx + 



-e 2 



dx 



+1 



(2m)' 



%E 



(A. 



(x+-ie) 2+1 J (x~-u)- 

This is the essential point where we assume that the dimension d is even. Otherwise the 
singularities appearing in these integrals become branch cuts and the residue theorem can 
not be applied. Hence, our final result is given by 



f2(E) 



d+2 

7T 2 



Cn 



(d+i) r(d-i)r(s±?) V2 



E 



(A.9) 



Next we turn to computing eq. (|A.4|) . Again we rely on the results of |27], where 
the form of the three-point function of the stress tensor was completely fixed using confor- 
mal invariance and energy conservation - in particular, see (3.15) in pTfl . 11 Inserting the 
appropriate indices and taking the limit xf — > oo, we find 



lim 



where 



d-2 



{T tj ix)T—( Xl )T i3 (0)) 



(_l)l-fd 2 -( d + 3 ) ( x -)2 t( x ) 



d+2 



d + 2 



d+6 



x + ie) 2 (x 1 — ie) 2 (x 2 ) 2 

(A.10) 



t(x) 
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^fix 3 ) 2 - x 2 ((x £ ) 2 + (x J ) 2 ) +4C 16(x J ) 2 (x j ) 2 - 2((x ? ) 2 + (x J ") 2 )x 2 + (x 2 ) 2 



+ [(d + 2){d- 2)A + {d- 2)B - 4dC] x 2 ((x J ) 2 + (x j ) 2 ) - 8(x 4 ") 2 (x i ) 2 
+ 8(d + 2){{d-2)A + B- AC) (x i ) 2 (x i ) 2 + 2(d + 4)(d - 2){B - 2C)(x i ) 2 (x i ) 2 . (A.ll) 

In the expression above, the ie prescription was t\ — > t\ — ie and t — )• t — 2ie, giving a larger 
negative imaginary part to time in operators standing to the left. Integrating over x^~ includes 
only one of the poles in eq. ( |A10 ), irrespective of whether the contour is closed in the upper 
or lower half of the xj~-plane. In either case, the integral yields 



dxi lim 
1 



d-2 



(Ii i (x)T__(x 1 )T ii (0)) 

r(d+l) 2-( ci+3 )i(a 



2vri 



rm 2 



2ie 



x d-l 



d+6 
(X 2 ) 2 



(A.12) 



11 However, we express our results in terms of the parameters A, B, C of jSfij. The corresponding parameters 
in Q are given by: a = A/8, b= (B — 2A)/8 and c = C/2. 
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Performing the remaining integrals over x a in eq. ( A. 4 ) as described in the previous case, the 
final result becomes 



h{E) 



E (E 



d+4 

2tt^ 



-5 ((d-2)(d + 2)A + 2dB-AdC) 



(A.13) 



Note that in this case, the integrals over the light-cone coordinates are 



dx + 



(x+ - 2ieY 



dx 



(x 



2ie 



.d+l 



(2tt«) 



i 1 fiE 
T{d+1) \ Y 



d+l 



(A.14) 



and so yield poles irrespective of whether d is even or odd. 

Combining the expressions in eqs. ( [All ) and flQjP then yields 

*V-2 h(E) d + l (d-2)(d + 2)A + 2dB -4dC 



(A.15) 



E f 2 (E) d (d-l)(d + 2)A-2B-4(d + l)C 

Note that producing this expression requires the use of an identity for T functions of the form: 

d-r 



n 2 2 - d T(d-l) 



(A.16) 



Comparing the previous expression with eq. ( |A.2[) , we find our first linear combination of ti 
and £4: 



1 



1 



-U 



d+l (d-2){d + 2)A + 2dB - idC 



(A.17) 



d-1" d 2 -!"* d (d-l)(d + 2)A-2B-4(d+l)C 

Following the same basic steps as above, one can also compute the flux with the vector channel 
polarization to produce the linear combination of ti and £4 in eq. (A. 3). The computations 
are somewhat more involved in this case and we only present the final result: 



d-S 

1 + — 7T*2 



1 , (d-2)(d + 2)A+ (3d-2)B -8dC 
2(d-l)" z d 2 ^l h ~ { + ' (d-l){d + 2)A-2B-4(d+l)C 

Combining eqs. ( |A.17| ) and ( |A.18[ ), we can solve for t% and t± separately, 

2{d + 1) (d - 2) (d + 2) (d + 1) A + 3 d 2 B - 4 d (2d + 1) C 



t-2 



d (d-l)(d + 2)A-2B-A(d+l)C 

(d + 1) (d + 2) {2d 2 -3d-3)A + 2d 2 (d + 2)B-4d{d+l)(d + 2)C 



(A.18) 



(A.19) 



d 



(d-l)(d + 2)A-2B-4(d + l)C 



Note that these expressions for general d reproduce the results presented in [18| and [17] in 
d = 4, 6 respectively. 
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B. Conformal tensor fields 



In sections 3J3 and [7|, we compare various results for the strongly coupled CFT's dual to GB 



gravity to those for massless free fields. As well as conformally coupled scalars and massless 
fermions, an (n-l)-form potential in d = In dimensions yields another free field theory which 
is also conformally invariant. The most familiar example is provided by d = 4, in which case 



d^xy/g^g""* ■ ■ ■ g^ F^ 2 ..., n F VlV2 ... Vn , (B.l) 



the potential is just an Abelian vector field. Our final results given in eq. ( B.15 ) are general 
expressions for A, B and C for any value of n. The analysis leading to these expression is 
quite involved and we only provide some of the salient steps in the following. Further, for 
simplicity, we work in a Euclidean-signature space. 

We begin with a free (n — l)-form gauge field A in d = 2n dimensions. As usual, the 
corresponding field strength is an n-form given by F = dA. With a general background 
metric, the action for this system is 

/= 2fr!) 

which is invariant under conformal rescalings, g^ v — > e 2 ^g M „. To simplify the subsequent 
calculations, we now fix the background metric to be simply flat space, i.e., g^ v = 5^ u . The 
corresponding energy momentum tensor is given by 

( n _l)l 'W-an^ 2(n!) v 

Because the theory is invariant under gauge transformations 5A = dA, in principle, 
we should introduce gauge-fixing terms to the action ( |B.1| ). However, correlation functions 
involving only gauge invariant operators will be unaffected by the choice of gauge-fixing and 
any ghost fields and so we do not consider the complete details here. With the natural 
generalization of Feynman gauge, the two-point function of the potential becomes 

r ( ft 1 ) i 

(A lll ... lln _ 1 (x 1 )A Ul ... Un _ 1 (x 2 )) = — 2 _ 1 Yl si g n (< 7 ) < W(ir-^™-i*> CT („-i) ( B - 3 ) 



T — F F a2 "' a " A F posv-an fB 21 



n-1 



where S n -i is a permutation group of n — 1 elements 1, 2, .., n — 1. Differentiating the above 
two-point function with respect to x^ n and x^" and subsequently antisymmetrizing the result 
with respect to u J s and z/s yields the following basic gauge invariant two-point function of 
the field strength 

where 

J Hi— tin ; "l—Un ( x ) = si S n (°") ^i^ CT (i)( x ) • • • ^nu^^ix) . (B.5) 



As in eq. (3.3), we have defined the orthogonal matrix 



V(x) = V - 2^ . (B.6) 
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After a lengthy calculation, we arrive at the connected two-point function of the energy 
momentum tensor 

n 2 T\2n — ll 1 

(T^(x)T^(O)) = ^ 2n J -^WW , (B.7) 



where T^^pix) is defined as in eq. (|3.2|). Of course, our result matches the expected form 
given in eq. ( |3.1| ) where 

C T = ^T[d-l] . (B.8) 

Proceeding further requires calculating the three-point function of the stress tensor. 
While such a computation is straightforward, it would be extremely tedious to carry out 
in general. The calculations are greatly simplified by going to the collinear frame [27|, for 
which the three points are chosen to lie on a straight line. In this case, the form of the 
three-point function is given by eqs. (4.20) and (4.21) of [27]. Since we wish to fix the values 
of A, B and C, it is enough to supplement eq. ( jB.8| ) with two extra relations between these 
constants. We choose to compute coefficients a and 7 defined in eq. (4.21) of j p7| . 

For simplicity, we exploit the rotational and translational symmetry of the problem to 
choose the following three points for the insertions of the energy momentum tensor: x\ = 
(x, 0), X2 = (y, 0), X3 = 0. As a result, 1^ becomes a constant orthogonal matrix given by 

I^v = <W - 2 SipSiv . (B.9) 

With some work, we find 



and hence, in the notation of 



d 2 r(d/2)r(d) 

a = 16 ^ • 

Similarly, we obtain 

(T^T^TniO)) = ^vrt L X 2n y 2n { ' _ y) 2n > (B.12) 

which yields 



d 3 T(d/2)T(d-l) 
7 " 3^W^ • (R13) 



According to eqs. (4.22-4.25) of |27f, 

(d-l)(d 3 -2d 2 -d + 4) (d-lf (d + l)(d-l)2(d-4) 

2d 2 Ad 2d 2 ' 

d'-d' 2d + 4 ^ + ( 1 _3) B _ (d+1)(d _ 1 ) c 

8 <i 8 2 <i 
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Hence combining these two equations with (|3,5| ),(B.8),( B.11|) ,( [B.13 ) and solving for A, B and 
C yields 



A 



d 3 T(d - l)r(d/2) 
8(d - 3)vr M /2 
B = (d-2)A , 

C = ——A . 



(B.15) 



These general results agree with those derived for d = 4 in [27] and for d = 6 in [17|. 



C. Useful expressions 

In order to write our equations as neatly as possible, we define the following expressions, 
which appear repeatedly in our analysis: 

/C = nV = 1-/ + A GB / 2 , 
^ = 1-2|§^A GB -2|^|A GB /(1-W), 
M = [1 - 2X GB f] 2 , 



V 



-m-i) 



1 ^-^ A 

1 " (£, _ 3) 



+ 



2(£>- !)(£>- 6) 
1+ (D^3) 



/- 



1 



[172 - 91L> + 1QD 2 - D 3 + 8(L> - !)(£> - 4)(2£> - 9)A GB ] A GB / 2 + 



2(D-3)(D-4) 

+ 2 ^ [87 - 44D + 5D 2 + 4(£) - !)(£> - 3)(£> - 4)A GB ] A GB / 3 - 



+ 



(£>-3) 2 

(D-5)p + 15) 3 4 8(D - 5) 4 5 
(D - 3) gb/ + (D — 3) gb/ ' 
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n = 



64/oo 



4/oo »' 
+ 4 
+ 



1 2 ^-^ A 
1 — — — a gb 



(D-3) 
2(D-1) 



2/ootx> 2 



(D-3) 



[(D-3)(D-4) 

A GB [4(£> - 1)A GB - 3D + 11] 



A GB [2(D - 1)A GB + D - 6] - 1 
/- 



+ 



+ 



2q 2 [4(D - 1)(2D - 7)A GB - 3D(D - 7) - 42] + 



(D-3)(D-A) 
+ 2(D - 4)/ 0O t» 2 A GB [4(D - 1)A GB - 7D + 31] 



A GB / 2 + 



+ 



(D-3)(D-A)[ 
+ q 2 [5D(D - 9) + 112] 



4A GB [q 2 (L> - 1)(2D - 7) - 2(D - 4)(£> - 5)/^] 
16(£> - 5) 



A 2 GB / 3 + 



(D_3)(D-4) 



[2q 2 (D-7)+ 



,. m , u m 2 A U 3 f 4 16p-5)(D-7)q 2 4 5 
+ 2(L» - 4)/ 00 w A GB J A GB j (£)-3)(£)-4) gbJ 



> = -^(D - 1) + ^ [3 + D + 10(D - 1)A GB ] / - ^(5L> + 11)A GB / 2 + i(£> + 7)A| B / 3 , 
r = _ (D-l) + 1 [5 _ D + g(jD _ i)AgJ / + 1 ( ^ _ 17)Agb/2 _ 1 (jD _ 9)a 2 b/3) 



V 



1 



(D - 1) + [D - 3 - 6(D - 1)A GB ] / - (D - 9)A GB / 2 + 2(L> - 5)A 2 B / 3 



2(D - 2) 

D. Coefficients for the sound channel 



Numerator of C. 



(i) 



2(D - 3)q 2 f 2 M 



2{D -1)(D- 2)X GB f O0 m 2 MMlC + q 2 AAP+ 

-(D- 3)q 2 7W 2 T + M\(D- 3) (q 2 P - (D - 2)f 00 vo 2 M 2 ) - <\ 2 MT 
-(D- 2)M 1/2 V \{D - l)q 2 /C - 2(D - 2)f 00 ro 2 M 1/2 ] [q 2 V - foo^M 

- q 2 / -2{D - lf{D - 2)A GB AW 1/2 /C 2 {f^M - q 2 V) + 

+ {D-1)(D- 3)MN 1/2 K \(D - 2)f (yo vo 2 M 2 + q 2 (MT -V) + 
+ (D — 2)q 2 (2M -Af)v]+ 2(D - 2) If^m 2 ((D - 3)M - N) NV+ 

+ q 2 [MV + (D-3)M 2 S]V 



+ 
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Denominator of c£Li =( D ~ Z)ufMN (fcc^N - q 2 V) 



2q 2 /((D-3)A4+A/) + 



(D - l)q 2 N 1/2 JC - 2{D - 2)f 00 xo 2 M 



Numerator of C. 



(2) _ 

sound 



(D-l) 



A/ 1 / 2 \(D - l)q 2 /C - 2{D - 2)f 00 v> 2 N 1 ' 2 



q 2 lC 



-{D-l) 2 uM 



-2q 2 V+ 



+ M {5(D - 2)f O0 xx> 2 M + 2q 2 T)] +{D-2) [M ((D - l)\ 2 uN - 8{D - 2)f 00 S) - 8(D - 2)/ 0O P] V 
+ A{D -!)(£>- 2)/ co q 2 A/ 1/2 /C 2 [V -M (T + 2{D - 2)V)] + 2{D -1)(D- 2)/ oc mx> 2 A/ 3/2 x 



+ 



+ 4q 2 / 2 



8 (D- l)(D-2) 2 / 00 A GB A^/C 2 x 



[-2q 2 7> + M[(D- 2)f O0 xx> 2 N - q 2 (25 + (D — 2)V)]] 

(D - 3)f OQ w 2 MM + q 2 ((D - 3)M + A/) vj + (D - l)q 2 WV 1/2 ((D - 3)7W + AT) x 
2q 2 A/7> — (D — 3)M 2 [(D - 2)/ 00 tx> 2 A/ + q 2 (2T - (D - 2)V)] + M [-(£> - 2) (2D - 5)/ 00 t» 2 AA 2 + 
+ 2(£> - 3)q 2 7> + q 2 A/ ((£> - 2) V - 27")] 1 + 2(D - 2)/ 0O /C [2q 2 A/ 2 7> - 4(L> - 3) 2 ql/W 3 T+ 
+ q 2 AW[2(D - l) 2 A GB mx> 2 A/ 3/2 + 6(£> - 3)V - 2Af (T + (D - 2){D - 3)V)] + 
+ (D - 3)A4 2 \2{D - l) 2 A GB q 2 m» 2 A/ 3/2 - 2(£> - 2){D - tyf^Af 2 + 4(D - 3)q 2 7>+ 



2q 2 A/ (3T + (£> - 2)(Z> - 3)V) 



+ 4q 2 / 



4(£> - 1) 2 (D - 2) 2 /^A GB n> 2 AlA/ 3/2 /C 3 + 

+ 2{D - 1)(D - 2)/ oc mx> 2 A/ 3 / 2 [-2q 2 A/7> + {D - Z)M 2 [(D - 2)f O0 x» 2 N + q 2 (T - S - (D - 2)V)] + 
+ A^ [(Z? - 2) 2 / 0O tx> 2 A/ 2 - 2(D - 3)q 2 P + q 2 A/ (T — S — (D — 2)V)]] + (D - 1)(D - 2)/ oc A/ 1/2 /C 2 x 
4q 2 A/7> - 2{D - 3)q 2 A4 2 (3T + 4(D - 2)V) + A/t [2(D - l) 2 A GB q 2 m» 2 A/ 3/2 + 
- 2(£> -2){D- 3)/ootx> 2 A/ 2 + 6(U - 3)q 2 7> - 8(D - 2) 2 / 00 A GB t» 2 A/ 1/2 V - 4q 2 A/(T + (£> - 2)V)] 



+ 



2A/7> [((£> - 1) Vn - 2(D - 2) 2 /^o> 2 ) N - 2(D - 2) 2 / oc q 2 V] + (D — 3)M' 



[AT[-3(D - 1) 2 (D - 2)/ oc qWA/ — 2 ((D — 1) 2 q 4 n - 2(D - 2) 2 /> 2 ) T] + 

+ (D — 2)q 2 ((£> - l) 2 q 2 uA/ - 8(D - 2)/ 00 <S) V] + MA/"[-4(D - 1) 2 (D - 2) 2 /^A GB mx> 4 A/ 3/2 + 

- D(D - 1) 2 {D - 2)/ oc qWA/ 2 + 2(D - 3)((U - l) 2 q 4 u - 6(U - 2) 2 /^tx> 2 )7' - 4(D - 2) 2 / oc q 2 5V+ 

+ A/[-2T {(D - 1) 2 q 4 « - 2(D - 2) 2 /^a, 2 ) + (U - 2){{D - 1) 2 q 4 u + 4(£> - 2) 2 (Z) - 3)/> 2 ) V]] 
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2q 2 f((D-3)M+AT) 
+ {D- l)q 2 M 1/2 IC - 2(D - 2)f 00 x» 2 N , 
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